1 -1 . A method for simulating fluid flow within a mold cavity, the method comprising 

2 the steps of: 

3 (a) providing a surface representation for a three-dimensional volume 

4 associated with a mold cavity; 

5 (b) separating the surface representation into at least a first portion and a 

6 second portion, the first portion of the surface representation being associated with at 

7 least one section of the volume having at least one of (i) a substantially invariant 

8 thickness and (ii) a gradually varying thickness along a length thereof; 

9 (c) discretizing a first portion of a solution domain bound on an exterior 

1 0 thereof by the first portion of the surface representation; 

1 1 (d) discretizing a second portion of the solution domain bound on an exterior 

12 thereof by the second portion of the surface representation; 

1 3 (e) defining a plurality of interface elements for the solution domain that 

14 connect at least part of the first portion of the solution domain to at least part of the 

1 5 second portion of the solution domain; 

16 (f) obtaining values of at least one process variable for the first portion of the 

17 solution domain using a first set of governing equations; and 

18 (g) obtaining values of the at least one process variable for the second portion 

19 of the solution domain using a second set of governing equations. 

1 2. The method according to claim 1, wherein step (b) is performed automatically. 

1 3. The method according to claim 2, wherein at least one of step (c), step (d), and 

2 step (e) is performed automatically. 
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4. The method according to claim 2, wherein at least two of step (c), step (d), and 
step (e) are performed automatically. 

5. The method according to claim 1, wherein the surface representation is a surface 
mesh. 

6. The method according to claim 1 , wherein the volume represents a molded object. 

7. The method according to claim 1, wherein the volume represents a mold cavity. 

8. The method according to claim 1, wherein the first set of governing equations in 
step (f), the second set of governing equations in step (g), and a set of interface element 
equations are solved simultaneously, subject to initial conditions and boundary 
conditions. 

9. The method of claim 8, wherein the interface element equations link a portion of 
the solution domain described by governing equations in two spatial dimensions to a 
portion of the solution domain described by governing equations in three spatial 
dimensions. 

1 0. The method according to claim 1 , wherein the at least one process variable is 
selected from the group consisting of temperature, pressure, fluid velocity, stress, and 
fluid flow front position. 

1 1 . The method according to claim 1, wherein there are at least two process variables 
selected from the group consisting of temperature, pressure, fluid velocity, stress, and 
fluid flow front position. 

12. The method according to claim 1, wherein there are at least three process 
variables including temperature, pressure, and fluid velocity. 
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13. The method according to claim 1, wherein the method simulates fluid injection in 
the three-dimensional volume. 

14. The method according to claim 13, wherein the method further comprises 
determining a location of at least one injection point. 

15. The method according to claim 1, wherein step (a) comprises providing the 
surface representation from CAD system output. 

16. The method according to claim 15, wherein the CAD system output is in 
stereolithography format or IGES format. 

17. The method according to claim 15, wherein the CAD system output defines a 
surface mesh comprising polygonal elements. 

1 8. The method according to claim 17, wherein the polygonal elements are triangular 
elements or quadrilateral elements. 

19. The method according to claim 15, wherein the CAD system output defines a 
three-dimensional mesh. 

20. The method according to claim 19, wherein the surface representation is provided 
from a lattice of polygons that bound the three-dimensional mesh. 

21 . The method according to claim 15, wherein step (a) comprises using the CAD 
system output as a preliminary mesh and remeshing the preliminary mesh to provide the 
surface representation. 

22. The method according to claim 1, wherein step (a) comprises providing a surface 
representation comprising a mesh of polygonal surface elements. 

23. The method according to claim 22, wherein step (b) comprises defining two or 
more subsurfaces, each subsurface comprising at least one of the surface elements. 
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24. The method according to claim 23, wherein step (b) comprises determining 
element properties and nodal properties for each of the surface elements. 

25. The method according to claim 24, wherein step (b) comprises using at least a 
subset of the element properties and nod^l properties to classify each of the two or more 
subsurfaces according to curvature. 

26. The method according to claim 23, wherein step (b) comprises defining at least 
one surface loop, each comprising a connected subset of edges of the surface 
representation. 

27. The method according to claim 23, wherein step (b) comprises remeshing at least 
a subset of the two or more subsurfaces using a bisection algorithm. 

28. The method according to claim 23, wherein step (b) comprises determining which 
of the two or more subsurfaces are matched subsurfaces. 

29. The method according to claim 28, wherein each pair of matched subsurfaces is 
separated by a definable thickness. 

30. The method according to claim 28, wherein the first portion of the surface 
representation comprises at least a subset of the matched subsurfaces. 

3 1 . The method according to claim 23, wherein step (b) comprises determining which 
of the two or more subsurfaces are unmatched subsurfaces. 

32. The method according to claim 3 1 , wherein the second portion of the surface 
representation comprises at least a subset of the unmatched subsurfaces. 

33. The method according to claim 23, wherein step (b) comprises determining which 
of the two or more subsurfaces are edge subsurfaces. 
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1 34. The method according to claim 28, wherein step (c) comprises projecting at least 

2 one of the surface elements from one subsurface in a substantially perpendicular direction 

3 onto a matched subsurface thereof, thereby defining paired surface elements. 

1 35. The method according to claim 34, wherein step (c) comprises converting the 

2 paired surface elements into wedge elements. 

1 36. The method according to claim 1, wherein step (c) comprises automatically 

2 discretizing the first portion of the solution domain. 

1 37. The method according to claim 1, wherein discretizing in step (c) comprises using 

2 the first portion of the surface representation to define the first portion of the solution 

3 domain. 

1 38. The method according to claim 1, wherein step (c) comprises discretizing the first 

2 portion of the solution domain using wedge elements. 

1 39. The method according to claim 38, wherein at least one of the wedge elements 

2 comprises at least one solution grid point along a thickness thereof. 

1 40. The method according to claim 38, wherein at least one of the wedge elements is 

2 a discretely layered element or a continuously layered element. 

1 41 . The method according to claim 1, wherein step (d) comprises automatically 

2 discretizing the second portion of the solution domain. 

1 42. The method according to claim 1, wherein step (c) comprises discretizing the first 

2 portion of the solution domain using hexahedral elements. 

1 43. The method according to claim 1, wherein step (c) comprises discretizing the first 

2 portion of the solution domain using shell elements. 
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1 44. The method according to claim 1 , wherein step (d) comprises discretizing the 

2 second portion of the solution domain using polyhedral elements. 

1 45. The method according to claim 44, wherein the polyhedral elements are 

2 tetrahedral elements or hexahedral elements. 

1 46. The method according to claim 1 , wherein step (e) comprises defining a set of line 

2 interface elements. 

1 47. The method according to claim 46, wherein each of the line interface elements is 

2 located along an interface of the first portion of the solution domain and the second 

3 portion of the solution domain. 

1 48. The method according to claim 46, wherein each of the line interface elements 

2 comprises at least two nodes of a wedge element of the first portion of the solution 

3 domain. 

1 49. The method according to claim 48, wherein each of the line interface elements 

2 further comprises at least one solution grid point between two of the at least two nodes. 

1 50. The method according to claim 1, wherein step (e) comprises defining a set of 

2 planar interface elements. 

1 51. The method according to claim 1 , wherein step (c) is initiated before step (e). 

1 52. The method according to claim 1 , wherein step (e) is initiated before step (d). 

1 53. The method according to claim 1 , wherein the first set of governing equations 

2 describes fluid flow in two spatial dimensions. 

1 54. The method according to claim 1 , wherein the first set of governing equations 

2 describes fluid flow in two spatial dimensions and time. 
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1 55. The method according to claim 1 , wherein the first set of governing equations 

2 describes fluid flow in one spatial dimension and time. 

1 56. The method according to claim 1, wherein step (f) comprises using a Hele-Shaw 

2 approximation. 

1 57. The method according to claim 1 , wherein step (g) comprises solving a Navier 

2 Stokes equation. 

1 58. The method according to claim 1 , wherein step (g) comprises solving a simplified 

2 Stokes equation. 

1 59. The method according to claim 1 , wherein the second set of governing equations 

2 comprises conservation of mass, conservation of momentum, and conservation of energy 

3 equations. 

1 60. The method according to claim 1, wherein at least one of step (f) and step (g) 

2 comprises using a meshless scheme. 

1 61 . The method according to claim 60, wherein the meshless scheme is a boundary 

2 element method, natural element method, or smooth particle hydrodynamics method. 

1 62. The method according to claim 1, further comprising the step of: 

2 (h) displaying the values of the at least one process variable directly on a 3D 

3 representation of the volume. 

1 63. The method according to claim 1, wherein step (g) comprises using a Mini 

2 element formulation. 

1 64. A method for simulating fluid flow within a mold cavity, the method comprising 

2 the steps of: 
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3 (a) providing a surface representation for a three-dimensional volume 

4 associated with a mold cavity; 

5 (b) automatically separating the surface representation into at least a first 

6 portion and a second portion; 

7 (c) defining a solution domain for the three-dimensional volume, where the 

8 solution domain comprises a first part corresponding to the first portion of the surface 

9 representation and a second part corresponding to the second portion of the surface 

10 representation; 

1 1 (d) solving for a process variable in the first part of the solution domain; and 

12 (e) solving for the process variable in the second part of the solution domain. 
1 65. The method according to claim 64, wherein the first portion of the surface 



2 representation in step (b) is associated with at least one section of the volume that has at 

3 least one of (i) a substantially invariant thickness and (ii) a gradually varying thickness 

4 along a length thereof. 

1 66. The method according to claim 64, wherein step (c) comprises automatically 

2 discretizing the first part and the second part of the solution domain. 

1 67. The method according to claim 64, further comprising the step of defining a 

2 plurality of interface elements that connect the first part of the solution domain to the 

3 second part of the solution domain. 

1 68. The method according to claim 64, wherein step (d) comprises using a first set of 

2 governing equations and step (e) comprises using a second set of governing equations. 

1 69. The method according to claim 68, wherein the first set of governing equations 

2 describes 2.5D flow and the second set of governing equations describes 3D flow. 
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1 70. A method for automatically defining a hybrid solution domain, the method 

2 comprising the steps of: 

3 (a) identifying a plurality of subsurfaces of a volume associated with a mold 

4 cavity using a representation of the surface of the volume; 

5 (b) matching one or more pairs of the plurality of subsurfaces to identify one 

6 or more matched pairs of subsurfaces and one or more unmatched subsurfaces; and 

7 (c) defining 

8 (i) a first portion of a hybrid solution domain bound at least in part by 

9 one or more of the matched pairs of subsurfaces and 

10 (ii) a second portion of the hybrid solution domain bound at least in 

1 1 part by one or more of the unmatched subsurfaces. 

1 71 . The method according to claim 70, wherein the volume represents a mold cavity. 

1 72. The method according to claim 71, further comprising using the hybrid solution 

2 domain to model a molding process. 

1 73. The method according to claim 70, wherein the representation of the surface of 

2 the volume comprises CAD system output. 

1 74. The method according to claim 70, wherein the first portion of the hybrid solution 

2 domain is amenable to 2.5D flow analysis, and the second portion of the hybrid solution 

3 domain is amendable to 3D flow analysis. 

1 75. The method according to claim 70, wherein step (b) comprises classifying each of 

2 the plurality of subsurfaces according to curvature. 

1 76. The method according to claim 70, wherein the matched pairs of subsurfaces each 

2 comprise two subsurfaces that are separated by a substantially constant thickness. 
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1 77. The method according to claim 70, wherein the volume represents a molded 

2 object. 

1 78. The method according to claim 77, further comprising using the hybrid solution 

2 domain in determining a structural property of the molded object. 

1 79. The method according to claim 78, wherein the structural property is warpage. 

1 80. An apparatus for simulating fluid flow within a mold cavity, the apparatus 

2 comprising: 

3 (a) a memory that stores code defining a set of instructions; and 

4 (b) a processor that executes said instructions thereby to 

5 (i) separate a surface representation of a three-dimensional volume 



6 associated with a mold cavity into at least a first portion and a second portion, the first 

i- 

7 portion of the surface representation being associated with at least one section of the 

8 volume having at least one of (i) a substantially invariant thickness and (ii) a gradually 

9 varying thickness along a length thereof; 



10 (ii) discretize a first portion of a solution domain bound on an exterior 

1 1 thereof by the first portion of the surface representation; 

12 (iii) discretize a second portion of the solution domain bound on an 

13 exterior thereof by the second portion of the surface representation; 

14 (iv) define a plurality of interface elements for the solution domain that 

15 connect at least part of the first portion of the solution domain to at least part of the 

16 second portion of the solution domain; 

17 (v) obtain values of at least one process variable for the first portion of 

18 the solution domain using a first set of governing equations; and 
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19 (vi) obtain values of the at least one process variable for the second 

20 portion of the solution domain using a second set of governing equations. 

1 81. An apparatus for defining a hybrid solution domain, the apparatus comprising: 

2 (a) a memory that stores code defining a set of instructions; and 

3 (b) a processor that executes said instructions thereby to 

4 (i) r identify a plurality of subsurfaces of a volume associated with a 

5 mold cavity using a representation of the surface of the volume; 

6 (ii) match one or more pairs of the plurality of subsurfaces to identify 

7 one or more matched pairs of subsurfaces and one or more unmatched subsurfaces; and 

8 (iii) define 

9 (A) a first portion of a hybrid solution domain bound at least in 

10 part by one or more of the matched pairs of subsurfaces and 

11 (B) a second portion of the hybrid solution domain bound at 

12 least in part by one or more of the unmatched subsurfaces. 

1 82. The method of claim 2, further comprising the step of re-characterizing a subset 

2 of the second portion of the solution domain as belonging to the first portion according to 

3 user input. 

1 83. The method of claim 2, further comprising the step of re-characterizing a subset 

2 of the first portion of the solution domain as belonging to the second portion according to 

3 user input. 

1 84. The method of claim 1 , wherein step (b) comprises separating the surface 

2 representation into a first portion, a second portion, and at least one additional portion. 
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